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We introduce a topological quantum number - coined dynamical topological order parameter 
(DTOP) - that is dynamically defined in the real-time evolution of a quantum many-body system 
and represented by a momentum space winding number of the Pancharatnam geometric phase. Our 
construction goes conceptually beyond the standard notion of topological invariants characterizing 
the wave-function of a system, which are constants of motion under coherent time evolution. In 
particular, we show that the DTOP can change its integer value at discrete times where so called 
dynamical quantum phase transitions occur, thus serving as a dynamical analog of an order parame¬ 
ter. Interestingly, studying quantum quenches in one-dimensional two-banded Bogoliubov de Gennes 
models, we find that the DTOP is capable of resolving if the topology of the system Hamiltonian has 
changed over the quench. Furthermore, we investigate the relation of the DTOP to the dynamics of 
the string order parameter that characterizes the topology of such systems in thermal equilibrium. 


I. INTRODUCTION 

In the theory of phase transitions mm, the phases sep¬ 
arated by a transition can be qualitatively distinguished 
in terms of order parameters which can either be local, 
following the conventional Ginzburg-Landau paradigm, 
or can be of topological nature. Regarding the latter, a 
complete classification for the ground states of gapped 
fermionic band structures in terms of global topologi¬ 
cal invariants has been achieved I5H5]. Out of thermal 
equilibrium, however, the current understanding of such 
systems is much less systematic. In particular, standard 
topological invariants [Si] associated with the instanta¬ 
neous wave function are constants of motion within uni¬ 
tary time evolution and can therefore not capture genuine 
non-equilibrium effects. In a non-equilibrium context, 
several recent studies |6H9| have discussed modifications 
to the bulk-boundary correspondence in non-equilibrium 
settings, indicating that the dynamical formation of edge 
states follows the instantaneous Hamiltonian and the sin¬ 
gle particle Green’s function, respectively. 

In this work, we construct a bulk topological quantum 
number - coined dynamical topological order parame¬ 
ter (DTOP) - which is dynamically defined: Going con¬ 
ceptually beyond the classification of topological phases 
in thermal equilibrium, it characterizes topological prop¬ 
erties of the real-time dynamics rather than of the in¬ 
stantaneous wave function or the instantaneous Hamil¬ 
tonian. In particular, we show in what sense the DTOP, 
somewhat analogous to order parameters at conventional 
phase transitions, distinguishes two phases separated by 
a nonequilibrium transition occurring in the coherent 
time evolution of a quantum system. Specifically, we find 
that the DTOP changes its value at a dynamical quan¬ 
tum phase transition (DQPT) [TU] which appears as a 
non-analytic behavior at critical times of the Loschmidt 
amplitude 

g{t) = {,p\m) = We-^^*\i’), ( 1 ) 

where \ip) denotes the initial state and H the Hamil¬ 


tonian governing the nonequilibrium quantum real-time 
evolution. Intriguing properties of DQPTs have been 
identified in a variety of different systems [TIJ1B5] . While 
Loschmidt amplitudes bear a formal similarity to equilib¬ 
rium partition functions uni, a dynamical analogue of an 
order parameter, that physically distinguishes the time 
intervals separated by a DQPT, has not yet been identi¬ 
fied. Here, we show how to construct such a dynamical 
analog for DQPTs by studying quantum quenches in two- 
banded Bogoliubov de Gennes models such as the Kitaev 
chain [21]. Notice that notions of dynamical transitions 
occuring out of equilibrum have also been introduced in 
different contexts [IS]. In the following, however, we will 
be referring to the definition made in Ref. m in terms 
of Loschmidt amplitudes (see Eq. 0 )- 

DQPTs occur whenever the time-evolved state |'0(t)) 
becomes orthogonal to the initial state vector |'0). This 
can be formally understood from the concept of parti¬ 
tion function zeros [2 [Ml [22], which occur as so called 
Fisher zeros in the context of Loschmidt amplitudes [lOj . 
Concerning the quest to identify dynamical order param¬ 
eters for DQPTs, this direct relation between DQPTs 
and wave function orthogonalities guides our intuition to¬ 
wards looking for an observable quantity that is smoothly 
defined for non-orthogonal state vectors and thus only 
allowed to behave discontinuously at critical times. The 
Pancharatnam geometrical phase (PGP) [Ml IM] is pre¬ 
cisely such a quantity. It was originally introduced |28j 
to define a relative phase for light beams with non- 
orthogonal polarization and has later been generalized 
to extend the notion of Berry’s geometric phase [Ml El] 
to general time evolution with non-orthogonal initial and 
final states, in particular allowing for non-adiabatic |32j 
and non-cyclic [M] dynamics. The geometric background 
of this construction is that a non-cyclic evolution can be 
augmented to a cyclic path in a unique way only if the 
two end states are non-orthogonal, namely by going back 
from the final to the initial state along a geodesic in pro¬ 
jective Hilbert space. 

The DTOP introduced in this work is a momentum- 
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FIG. 1. (color online) Left panel: Color plot of the Pancharatnam geometric phase see Eq. (j^, for chemical-potential 

quenches fj, = 0 ^ = 3 in the Kitaev chain as a function of lattice momentum and time. Time is measured in units of 

the critical time tc where the first dynamical quantum phase transition (DQPT) occurs. The critical momentum fcc at which 
non-analyticities occur is marked with a black dotted line. Right panel: Rate function g{t) = —Af“^Relog[|(J(t)p] of the 
Loschmidt echo C{t) = see Eq. Q, and the dynamical topological order parameter voit), see Eq. (|^, as a function 

of time. The real-time non-analyticities in g{t), occuring at odd multiples of tc (red dashed lines), define the DQPTs. unit) 
changes its value only at the DQPTs thus serving as a dynamical order parameter. 


space winding number of the PGP which serves as a 
dynamical analog of a topological order parameter in 
two-banded Bogoliubov de Gennes models undergo¬ 
ing a DQPT after a quantum quench, i.e., a sudden 
change in the band structure parameters. We show 
that the integer-valued DTOP can change its value 
only at DQPTs and, moreover, that this allows us to 
dynamically resolve how the topology of the underlying 
Hamiltonian has changed during the quench. Our con¬ 
struction relies on the presence of a so-called particle hole 
symmetry (PHS), i.e., a spectral constraint imposed by 
an antiunitary operator C, = 1 which anti-commutes 
with the system Hamiltonian. In Bogoliubov de Gennes 
models, such a constraint is naturally imposed by the 
fermionic algebra to the Nambu spinor representation of 
the Hamiltonian [33]. We illustrate our construction by 
studying quenches in several models including, e.g. the 
Kitaev chain [33] . 


Outline. The remainder of this article is organized as 
follows. Section El is devoted to the definition of our 
model system as well as the general analysis of its quench 
dynamics, in particular the occurrence of DQPTs. In 
Section [TTll we construct the DTOP which is at the heart 
of this study and discuss its basic properties. Thereafter, 
in Section [IV[ we present benchmark simulations, show¬ 
ing how the DTOP characterizes the dynamics of several 
ID Bogoliubov de Gennes models. The relation to the 
familiar string order parameter which hallmarks differ¬ 
ent phases in equilibrium is investigated in Section |V| 
Finally, in Section VI we outline how the physics dis¬ 


cussed in this manuscript, in particular the DTOP, could 
be experimentally observed, and put our present work 
into a broader context by discussing its relation comple¬ 
mentary approaches. 


II. MODEL AND DYNAMICS 

A. Underlying Hamiltonian 

The dynamical properties we are concerned with here 
are generated by gapped free fermionic two-banded Bo- 
goliubov de Gennes models in ID without requiring fur¬ 
ther symmetries, i.e., in symmetry class D |33j . We de¬ 
note the Nambu pseudo-spin by t and choose the conven¬ 
tion C — tiK for the PHS operation, where K denotes 
the complex conjugation. Assuming a unit lattice con¬ 
stant, the first Brillouin zone is the circle resulting from 
the interval [— 7r,7r] by identihcation of its end points. 
The Bloch Hamiltonian is then of the form 

3 

H{k)=d{k)-T = '^d^{k)Tj ( 2 ) 

i=i 

and satisfies the spectral PHS constraint 

T,H{k)n = -H*{-k). (3) 

As a consequence of Eq. ([^, d^{k) and d'^{k) must be odd 
functions of the lattice momentum k, while d^{k) must 
be even. Eq. (|^ is local in momentum at the two real 
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lattice momenta fcR = 0, tt which satisfy k = —k (mod27r) 
and where both d} and (P need to vanish such that 

H{kK) = PikR)T 3 , kR e {0,7r} . (4) 

With the unit vector d(k) = d{k)/\d{k)\, the geometrical 
interpretation of Eq. (|4j) is that d is pinned to the poles 
of the Bloch sphere at the real momenta which plays a 
crucial role in the following. There are two topologi¬ 
cally inequivalent classes of such Hamiltonians [23] , here 
simply distinguished by the sign of d^{0)d^{7r) which be¬ 
comes negative for the non-trivial topological phase. 


B. Quench dynamics and DQPTs 

We study nonequilibrium quantum real-time evolution 
and DQPTs induced by a quantum quench. The system 
is prepared in the ground state lijj) of an initial Hamilto¬ 
nian Hi(k) = diik) ■ T. At time t = 0, a parameter will 
be switched suddenly within the set of models in Eq. Q 
resulting in a sudden change di{k) >-)■ df{k). We assume 
that the system initially occupies the lower Bloch band 
of Hi. The associated lower Bloch states are denoted 
by such that the initial state IV') is a Slater de¬ 

terminant of all lower band Bloch states. Since lattice 
translation invariance is maintained at all times, the dy¬ 
namics of the system can be considered separately for 
every lattice momentum k. Explicitly, we get 

IV'feW) = (5) 

where ±e{ = ±\df{k)\ denotes the energy eigenvalues of 
Hfik), its Bloch states, and gk = = 

with Ifffep = i(l -h d^k) ■ df{k)), \ek\^ = 
|(1 — di{k) ■ df{k)) are expansion coefhcients of the ini¬ 
tial lower Bloch state in the new Bloch states after the 
quench. For the geometric interpretation of our later re¬ 
sults, it is helpful to consider the vector dpk) as a refer¬ 
ence direction, say pointing to the south pole of a Bloch 
sphere defined at every momentum, and to consider the 
direction of df{k) relative to this reference. We refer to 
this construction as the relative Bloch sphere in the fol¬ 
lowing. 

Recently, DQPTs in topological systems satisfying 
Eq. as well as in related spin chains have been identi¬ 
fied T3[T7||2T]. DQPTs are caused by Fisher zeros [TO] 
where for a momentum k^ the overlap = 0 

vanishes at a time tc- Here, this can only happen if 

IgfceP = |e/ccP and n € N. (6) 

24. 

Fisher zeros and DQPTs hence occur at momenta where 
the initial lower Bloch state is an equal weight superpo¬ 
sition of the final Bloch states, i.e. at di{kc) ■ df{kc) = 0 
marking the equator of the relative Bloch sphere. 


whereas the critical time is determined by the spectrum 

f 

ej. of the final Hamiltonian. 


III. CONSTRUCTION OF A DYNAMICAL 
TOPOLOGICAL QUANTUM NUMBER 

A. Pancharatnam geometric phase 

In order to define the PGP at lattice momentum 
fc, let us decompose the Loschmidt amplitude G{t) = 
l\k>oSk{t) with 

5k{t) = {ul-\Mt))=rkit)P'^'‘^*\ (7) 

and rk(t), 4’k{t) its polar coordinates. The phase (jjkit) 
contains a purely geometric and gauge-invariant compo¬ 
nent 


= Mt) - , ( 8 ) 

obtained by subtracting the dynamical phase = 

-/o ds{ipkis)\Hf\'ipk{s)) = efctdfffcp - |efcp). (I)k is the 
aforementioned PGP [25] that will be the central build¬ 
ing block for the DTOP introduced in this work. We 
stress that this definition of the PGP becomes singular 
at Fisher zeros as the total phase in Efi- 0 is 

ill-defined at critical times. 


B. Definition of the DTOP 

Eq. Q implies that either |efcp = 0 and |(?fcp = 1 or 
vice versa at the real momenta /cr = 0, tt. From Eq. (§, 
we directly conclude that (l>ks{t) = i.e., that the 

PGP is pinned to zero at these special momenta. Thus, 
as far as the PGP is concerned, the interval [0, tt] between 
the real momenta can be endowed with the topology of 
the unit circle by identifying its end points. We refer 
to this periodic structure as the effective Brillouin zone 
(EBZ). We are now ready to define a DTOP in terms of 
the PGP as 


i^oit) 


^ rMw. 

27r /o dk 


(9) 


is the integer-quantized winding number of the 
PGP over the EBZ and is smoothly defined as a function 
of time in the absence of Fisher zeros. More formally, 
is a topological invariant distinguishing homotopi- 
cally inequivalent mappings EBZ —>■ U{1), k 1 —>■ 
from the unit circle to itself. The definition of voit) 
in Eq. (|^ and its subsequent further interpretation are 
the main results of our present work. 

Since DQPTs can only occur at points in time where 
Fisher zeros are present, I'nit) rnust be constant in time 
intervals between DQPTs as it cannot smoothly change 
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its integer value. But does vd change its value at ev¬ 
ery DQPT? We answer this question in the affirmative 
implying that vd (t) can serve as an order parameter for 
the studied DQPTs. We find that, quite remarkably, the 
change in the DTOP Avu{tc) in the vicinity of a critical 
time tc can be directly related to the sign of the slope 
Ska ~ ('9fe|6fcP)|fe the critical momentum as 

Avoitc) = lim [vnitc + t) - voitc - t)] = sgn(sfcj 

T —^0 

( 10 ) 


which loosely resembles an index theorem. This result 
affords an intuitive geometric interpretation: As pointed 
out before, critical momenta are located on the equator 
of the relative Bloch sphere. Avo{tc) is then directly 
related to whether df{k) traverses the equator of the rel¬ 
ative Bloch sphere from the northern to the southern 
hemisphere (sgn (s^^) = —1) or from the southern to the 
northern hemisphere (sgn(sfc^) = 1) at the critical mo¬ 
mentum. 

To establish Eq. (10), we first identify a fundamen¬ 
tal dynamical symmetry of Qk{t) (see Eq. 0 ) at crit¬ 


ical momenta k^- From Eqs. 

Gkait) G Furthermore, the dynamical phase is zero 


5][6), we conclude that 


= sgn (^cos{elj 


due to I P = I P such that e*'^''c = : 

i.e., the PGP is pinned to the real values 0, tt at the criti¬ 
cal momenta at all times. When passing through critical 
times, marked by cos{eljc) = 0 , the sign of cos(e{^t) 
changes and the PGP jumps by tt. Expanding dk(l>k{t) 
around kc and tc to leading order, it is straightforward 
to prove Eq. (10). We note that the connection between 
Fisher zeros and tt jumps of the PGP is generally valid 
beyond the scope of the present DTOP: When a complex 
function (the Loschmidt amplitude) goes through zero as 
function of a real parameter (time), its phase jumps by tt. 
Since the dynamical phase is always continuous in time, 
this jump can only occur in the PGP. 


IV. BENCHMARKS SIMULATIONS 

We further investigate and illustrate the DTOP voit) 
with two benchmark examples. First, we study the Ki- 
taev chain [M]. a toy model for a proximity induced p- 
wave smerconductor. The model Hamiltonian is of the 
form ([2) with d}-{k) = 0, c?^(fc) = sin(fc), (Pik) = 
pL — cos(fc), where physically the cos{k) term represents 
the kinetic energy and the sin(fc) term represents the p- 
wave pairing. Here, we simulate a quench in the chem¬ 
ical potential from /i = 0 in the initial Hamiltonian 
Hi to /X = 3 in the final Hamiltonian Hf. At the 
real momentum k = 0 the (P component switches sign 
over the quench while at fc = tt it does not, mean¬ 
ing that the topological phase of the Hamiltonian [24] 
changes during the quench. Hence, |eoP = 1 while 
|e 7 rP = 0. Due to continuity, there must be a crit¬ 
ical momentum kc in the interior of the EBZ where 


IcfecP = IfffecP = 5 in agreement with Ref. jH]. Accord¬ 
ing to Eq. (1^, this implies that DQPTs will occur at 
times = (2n — 1)tt/{ 2eka) — — l)tc; n G N. More 

specifically we find, kc = arccos(l/3), tc = 7r/(4-\/2) 
and Ska = “Ij i-®-) crosses the equator of the 

relative Bloch sphere in southern direction, correspond¬ 
ing to a change A{tc,n) = — 1 in the DTOP (see Eq. 
([I0| )) at the critical times tc,n = (2n — l)tc- In Fig. 
[ipwe show a color-plot of the PGP as a function of k 
and t from which the critical times and the associated 
changes in the phase winding number vd representing 
our DTOP become visually clear (left panel). Moreover, 
Fig.g displays the time dependence of both the DTOP 
VD{t) and the rate function g{t) = —N~^ log(|C/(t)p) = 

-7r“^Re (^log/J" d/c IfffeP-H ^ which plays 

the role of a thermodynamic potential here and whose 
points of nonanalytic behavior define the DQPTs m 
(right panel). The DTOP indeed changes its value at 
every DQPT and uniquely characterizes the dynamical 
phase in between two DQPTs. 

As a second benchmark, we consider a quench from 
di{k) = (0,0,1) to df{k) = (0,sin(A:),l -I- cos(2A:) -I- 
Acos(fc)). This model is similar to the Kitaev chain 
studied before, but also includes a next-to-nearest neigh¬ 
bor hopping. For 0 < A < 2, this quench does not 
change the topological phase of the Hamiltonian since 
di{k]g^) = df{ks)- Still, from Eq. ^ we find two crit¬ 
ical momenta ki^'^ = 7r/2, k!?'^ = arccos(—A/2). At 
ki^\ df{k) enters the northern hemisphere of the relative 
Bloch sphere and returns to the southern hemisphere at 
k^\ i.e., di{k)-df{k) < 0 for < k < k^\ As a conse¬ 
quence, sgn(s^(i)) = —sgn(s^( 2 )) = 1 and, from Eq. (10), 
we see that the change in the DTOP is opposite at the 
two critical momenta. In Fig. [^ we show a color plot of 
the PGP for this quench (left panel) as well as the time 
dependence of the rate function g{t) and the DTOP voit) 
(right panel). The DTOP changes at every DQPT, how¬ 
ever, due to the competing A{ti%) = (—1)*+^, z = 1,2, 
its behavior is not monotonous as opposed to the first 
quench example. 


V. RELATION TO STRING ORDER 
PARAMETER 

It is natural to ask how the DTOP defined in Eq. 
is connected to the underlying equilibrium topology of 
the system. In this context, we show in the following 
numerical evidence that the dynamics of the DTOP can 
be linked to the decay of string order parameters. In 
thermal equilibrium, the string order parameter directly 
reflects the bulk topological properties. Focusing on the 
continuation of such bulk properties to non-equilibrium 
systems, our present study complements previous work 
reporting deviations from the bulk-boundary correspon¬ 
dence out of thermal equilibrium ISHH]- There, surface 
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FIG. 2. Left panel: Color plot of the PGP as a function of lattice momentum and time for A = 1.3. Time is measured 

in units of the first critical time = vr/2. The critical momenta and at which DQPTs appear as discontinuities in 
the PGP is marked with black dotted lines. Right panel: Rate function g{t) and DTOP voit) as a function of time for A = 1.3. 
The DQPTs occur at odd multiples of the critical times and 


states have been shown to develop dynamically even if 
the conventional bulk topological quantum numbers re¬ 
main trivial [5]. 

Rather than exhibiting local order, topological phases 
are typically characterized by nonlocal properties. The 
topological phase in the Kitaev chain [53] (see also first 
benchmark in Section IV), is associated with a nonva¬ 


nishing expectation value of the so called string order 
parameter 


Ol^ = {ci + {cm + cl). (11) 


Here, ci with I = 1,..., N denotes the fermionic an¬ 
nihilation operator in its real-space representation. In 
Fig. we show numerically obtained dynamics 0{t) = 
lim|;_,,„l_,,oo(0;m(i)) for the string order parameter in 
the Kitaev chain for the same set of parameters as in 
Fig. HI but for a slightly more general model including 
also a next-to-nearest neighbor hopping with strength j, 
i.e., d^{k) = gL — cos{k) — jcos(2fc). We have calculated 
the dynamics of O using the mapping onto Pfafhans [Mj 
which can be evaluated numerically very efficiently. For 
j = 0 where the Kitaev chain can be mapped onto a 
transverse-field Ising model, the string order parameter 
maps to the spin-spin correlation function of the Ising 
order parameter whose dynamics has been studied pre¬ 
viously [TUI 155] . 

As one can see in Fig. starting initially in the topo¬ 
logical phase and quenching into the trivial phase, the 
string order parameter exhibits a damped oscillatory de¬ 
cay. The decay of the string order parameter is, of course, 
compatible with the topologically trivial structure of the 
final Hamiltonian. In the context of the DTOP, it is, 
however, of particular importance how the string or- 



FIG. 3. Decay of the string order parameter for the same 
chemical potential quench as in Fig. [T] but for a model also 
including a next-to-nearest neighbor hopping with amplitude 
j. The time axis is rescaled by t* where DQPTs occur and 
the DTOP changes its quantized value. The inset shows a 
comparison of the low-energy time scale tgap(i) = 
with A(j) the gap, and tc{j). The simulations are done for a 
lattice system with up to A = 100 sites, where we find that 
the data is converged with respect to system size. 


der decays. Specifically, the time scale for the oscilla¬ 
tions on top of the exponential envelope coincides with 
tc{j) setting the critical times where DQPTs occur and 
the DTOP changes dynamically. To show this, we have 
rescaled the time axis by tc{j). After an initial tran¬ 
sient regime where the time scale for the oscillations does 
not perfectly match tc{j) yet, in the asymptotic decay 
regime the agreement becomes very good. To further 
strengthen this suggested connection, we have also in- 
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eluded an inset where we compare the emergent nonequi¬ 
librium time scale tc(j) with the conventional time scale 
tgapij) = A(j), i.e., the gap of the j-dependent 

the final Hamiltonian. tc{j) we have extracted by solv¬ 
ing Eq. (|^. While tgap(j) increases with increasing j, we 
find that tc{j) decreases, thereby excluding the possibil¬ 
ity of an accidental similarity of the two scales. This gives 
further evidence that the time scale associated with the 
string order parameter is not set by the gap but rather 
by the emergent nonequilibrium time scale tc{j) and thus 
is directly connected to the dynamics of the DTOP. 

VI. CONCLUDING DISCUSSION 

The DTOP vd defined here (see Eq. Q) qualitatively 
distinguishes periods of time-evolution that are separated 
by DQPTs. As a truly dynamical quantity, the DTOP is 
fundamentally different from the conventional topological 
invariants that classify ground states of gapped Hamil¬ 
tonians HU. In particular, the conventional topological 
invariant for the time-dependent state \tjj(t)) is a, constant 
of motion in our present non-equilibrium setting [51l^ll6j. 
However, there is an interesting interplay between equi¬ 
librium invariants and the occurrence of DQPTs. In this 
context, it has been shown |2Ij that a quench between 
topologically inequivalent Hamiltonians implies the pres¬ 
ence of DQPTs, in agreement with the situation in our 
first benchmark example. However, DQPTs can also hap¬ 
pen if the initial and the final Hamiltonian are topolog¬ 
ically equivalent as in our second benchmark example. 
Remarkably, the structure of the DTOP is capable of re¬ 
solving these different scenarios. From Eq. ( [l0| ), we see 
that the DTOP behaves qualitatively different in these 
two cases (see comparison of right panels of Fig. [^and 
Fig. §: If the Hamiltonian topology changes, there is an 
odd number of critical momenta giving rise to a change 
of the DTOP after one DQPT associated with each of 
the critical momenta. In contrast, for quenches between 
equivalent Hamiltonians, the sum over the changes of vd 
for all critical momenta is zero. 

The construction of the DTOP relies on the two- 
banded character of the models we consider. For larger 
unit cells, pertaining in particular to the modelling of dis¬ 
ordered systems, the change of the PGP between the real 
momenta is not necessarily quantized because the spec¬ 
tral PHS constraint in Eq. ^ does not enforce Eq. Q. 
However, in many physical situations considering an ef¬ 
fective model with two bands is a physically well justified 
approximation and the DTOP defined here is hence also 
expected to emerge in more complex systems. Further¬ 


more, the general observation that Fisher zeros in real 
time lead to tt phase jumps in the PGP is valid beyond 
the specific construction of the DTOP. This may serve 
as a starting point for the generalization of our present 
construction to higher spatial dimensions in future work. 

The dynamical order parameter identified in this work 
is of topological nature. A natural question is whether 
also local dynamical order parameters can exist. Al¬ 
though a general answer is beyond the scope of this work, 
we can provide some intuitive insights. A dynamical ana¬ 
log of a diverging correlation length, giving rise to nonan- 
alytic behavior of local observables in continuous equilib¬ 
rium phase transitions, cannot develop dynamically in a 
finite period of time. This is due to fundamental causal¬ 
ity constraints such as Lieb-Robinson bounds. Hence, 
we do not expect direct analogs of such phenomena at 
DQPTs. 

We conclude by summarizing some recent experimental 
progress on the building blocks of a setup where DQPTs 
and the DTOP could be observed. Models similar to 
the Kitaev chain [24] . see our benchmark examples, can 
be realized both in solid state systems and po¬ 

tentially with cold atoms in optical lattices [121 El] • In¬ 
stead of actually realizing a superfluid system described 
by a Bogoliubov de Gennes equation, experimental stud¬ 
ies may resort to insulating band structures similar to 
the Su, Schrieffer, and Heeger model [IS] which also ex¬ 
hibits a formal PHS and for which the DTOP is analo¬ 
gously defined. While inducing nonequilibrium dynamics 
in solid state systems is challenging, with ultracold atoms 
quantum quenches have already been studied experimen¬ 
tally (Ml EH- Moreover, momentum-resolved phase dif¬ 
ferences of Bloch functions, needed for the DTOP, have 
recently been measured in terms of Berry phases [151150] . 
In particular, the experimental techniques employed in 
Ref. [50] can be directly employed to extract the PGP 
and to reconstruct the DTOP defined in the present 
work. Paving the way towards the general observation 
of DQPTs, a measurement scheme for Loschmidt echos 
= |^(t)P has been introduced [5T1 [52] . 
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